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THE GAMMA FUNCTION AND APPLICATIONS 
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Abstract. In the article we present necessary and sufficient conditions for 
a function involving the logarithm of the gamma function to be completely 
monotonic and apply these results to bound the gamma function T(x), the 
n-th harmonic number 5Zfe=i x> an£ I the factorial n\. 



1. Introduction 

We recall from [18, Chapter XIII], [29, Chapter 1] and [31, Chapter IV] that a 
function / is said to be completely monotonic on an interval I if / has derivatives 
of all orders on / and 

< (-l) n f (n) (x) < oo (1.1) 

for x € I and n > 0. 

We also recall from [1, p. 254, 6.1.1] that the classical Euler gamma function 
T(x) may be defined by 



T(x) 



t*-Vdt, 



x>0. (1.2) 



_ r'Qr) 

r(x) 



is called the psi or 



The logarithmic derivative of T(x), denoted by ip(x) 

di-gamma function, and the derivatives tp^(x) for i £ N are respectively called 
the polygamma functions. They are a series of important special functions and 
have much extensive applications in many branches such as statistics, probabil- 
ity, number theory, theory of 0-1 matrices, graph theory, combinatorics, physics, 
engineering, and other mathematical sciences. 

Using Hcrmitc-Hadamard's inequality (see [26, 27, 28]), the double inequality 



1 

X -~2 



\n\x-- 



< x 



In x 



In V2vr 
- 1 



1 



24(x- r 
lnV^ 



< lnT(x) 



24(^r 2 +x + l/2 - 1/2) 



for x > 1 was obtained in [4, p. 236, Theorem 1] 
In [30, p. 1774, Theorem 2.3], the function 



H(x) = lnr(a:+l) 



24(1 + 1/2) 



(1.3) 



(1.4) 
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was proved to be completely monotonic on (0, oo). From this it was deduced in [30, 
p. 1775, Corollary 2.4] that the double inequality 



, 1 ,r > \x + l/2 
X + i / Z \ „-l/24(x+l/2) 



e 



</3 ' g-i/M^+i/a; (L5) 



holds for .t > 0, where a = V2n = 2.50 • • • and (3 = V2e 7 / 12 = 2.53 • • • are the 
best possible constants. 

We observe that, by taking the natural exponentials on all sides of (1.3) and 
replacing x by x + 1, the inequality (1.3) may be rewritten as 



< v ^( r £±VHy^ , %-V24(Vx='+3x+5/2-l/2) ) x > Q (L6) 



s+1/2 



Hence, it is clear that the left hand side inequality in (1.5) is stronger than the 
corresponding one in (1.3) or, equivalently, (1.6). But, when x > 1, the right hand 
side inequality in (1.5) is weaker than the corresponding one in (1.3) or, equivalently, 
(1.6). 

For A > and x e (0, oo), let 
H X{X ) = In V(x + 1) - (, + I) ln(* + I) + x + \ - In yfa + . (1.7) 

The first aim of this paper is to find necessary and sufficient conditions for the 
functions ±H\(x) to be completely monotonic on (0,oo). The second aim is to 
apply the complete monotonicity of ±H\(x) to establish inequalities for bounding 
the gamma function, the n-th harmonic number Y^k=i an d the factorial n\. 

2. Necessary and sufficient conditions 

Our main results are necessary and sufficient conditions for the functions ±H\(x) 
to be completely monotonic on (0, oo), which can be stated in the following theorem. 

Theorem 1. For x e (0, oo) and A > 7 

(1) if and only ifO < A < |, the function H\{x) is completely monotonic; 

(2) if 

"24/1 t 



A > — inf { - In . . 

~ t£(o,oo)\ t [t 2 Ve*/ 2 e* - 1 

the function —H\(x) is completely monotonic; 



(2.1) 



(3) if A > | 7 a special case of the inequality (2.1), the function —H\(x) is 
completely monotonic on (0, oo). 

Proof of the necessary condition in (1) of Theorem 1. If the function H\(x) is com- 
pletely monotonic on (0, oo) for A > 0, then H\(x) > 0, which may be rearranged 
as 
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where 



f(x) = lnT(x + l)-[x+^)ln[x+^-)+x+^-- ^ln(27r), x > 0. (2.2) 



Since 

as z — > oo in |arg2:| < 7r, see [1, p. 257, 6.1.41], we have 

^~H)'"^-^i2^ +o (>) < 2 - 4 » 

as x — » oo. As a result, we have 

1 1 fei(x) 

24/(x) ~ ~2 ' ^) 

as x — > oo, where 



and 



hAx) = 12x(2x 2 + 3x + 1) In - \ \- - I2x 2 - 9x + 1 + 0(1) -> ^ 
v ' x + 1/2 2 



2 , , • .'• • I .. . . V I \ 1 



h 2 (x) = 6(2x 2 + 3x + 1) In ^— -6x-5 + 0. , 

v 7 x + 1/2 \x J 2 



as x — > oo. So 

- x ~2Tf(x)^ l r x ^°°- 

This means that the necessary condition for the function H\(x) to be completely 
monotonic on (0, oo) is A < \. 

First proof of the sufficient condition in (1) of Theorem 1. When A < \, the 
function H\(x) is 

H x (x) = H 1/2 (x) + - 2 4(x+l/2) 

i r 1 ' 2 i 

Since the function Hi/2(x) = H{x) is completely monotonic on (0, oo), see [30, 
p. 1774, Theorem 2.3], and the integral j^ 2 ^jjp dt is clearly completely mono- 
tonic on (— A, oo ) D (0,oo), also since the product and sum of finite completely 
monotonic functions are completely monotonic on their common domain, it follows 
immediately that the function H\(x) is completely monotonic on (0, oo) if A < 1. 

Second proof of the sufficient condition in (1) of Theorem 1. The famous Binet's 
first formula of lnT(x) for x > is given by 

lnT(x) = (x- ^ lnx-x + lnv 7 ^ +9(x), (2.5) 

where 

9W =i, (—-- t + 2)— dt < M » 
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for x > is called the remainder of Binet's first formula for the logarithm of the 
gamma function, see [17, p. 11] or [25, p. 462]. The formulas 



T(z) = k z 



t z - x e- kt dt 



and 



a 



du 



(2.7) 



(2.8) 



for Rez > 0, Refc > 0, a > and b > can be found in [1, p. 255, 6.1.1 and p. 230, 
5.1.32]. Utilizing these formulas yields 



a; + 1/2 2 
'(*) 

oo/ e -t/2 1 



t 



i 



■<?(*), 



dt, 



H' x (x) = f'(x) 



24(x + A) 2 
e -t/2 ! 



-XI 



(2.9) 
(2.10) 

(2.11) 



" x *di. 



t e* - 1 24 
From (1.7), (2.2) and (2.9), it is easy to see that 

lim H\(x) = lim f(x) + lim — 1 = 0. 

x— 7-oc 2:^00 >oo Z4(X + A ) 

Therefore, in order to prove the complete monotonicity of =pH\(x), it suffices to 
show ±H' x (x) is completely monotonic on (0, oo). For this, it is sufficient to have 



(2.12) 



-t/2 



1 



te 



-At 



24 



< 







for alii € (0, oo), which is equivalent to 



A | --in 



24 fe^' 2 



t V t 



(2.13) 



for all t e (0,oo). 
We claim that 



In 



24 /e'*/ 2 



1 



t V t 



1 



> 



1 



for t G (0, oo). In fact, this inequality can be reduced to 

24 /e-'/ 2 1 



t V i 



1 



(2.14) 



equivalcntly, 

OO 

(t 2 - 24)e* + 24ie*/ 2 - t 2 + 24 = { [fc(fc - 1) - 24]2 fe + 48/c} ^ 



fc=5 
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n 00 4-k 

= ^ 5 +E{[ fc ( fc - 1 )- 24 ] 2fe+48fc } 



fc=6 



7 00 t k 

^ Y/ + E^( fc - 1) - 24](1 + fc) + 48k} m 



fc=6 



> ^t 5 + f> 3 + 23fc - 24 N ' 



24 
> 0. 



fc=6 



fc!2 fc 



Thus, when < A < |, the function H\{x) is completely monotonic on (0,oo). 

Proof of (2) in Theorem 1. This follows from the inequality with the sign > 
in (2.13). 



Proof of (3) in Theorem 1. Suppose that 

24 /e-*/ 2 



In 



for t e (0,oo). Then 



t V t e* - 1 

24 fe-t' 2 1 

> e 



< A 



(2.15) 



-At 



i V * e* - 1 , 
which can be rewritten as 

24e At (e t -te*/ 2 -l) > iV^e* - 1). 

Expanding at t = the functions on both sides of the above inequality into power 
series yields 



24 E 

fe=3 



(A+l) fe -A fe -fc( A+ 



fc-i 



f fe 00 

' ^E 



fc! 



fe=3 



fc-2 / ^ \ fc-2' 



(fc-2)! 



(2.16) 



Let hk;\{u) = (A + u) k 1 . Then, by the left hand side inequality in the double 
integral inequality 



(b - a) 2 1 
-to < 



24 



6 — a 



° + ^<t^M, (2.7, 



where / : [a, b] — > M be a twice differentiablc mapping and m < f"(t) < M for all 
t e (a, 6), see [5, 6] and [26, p. 236, Theorem A], we obtain 



(A + 1) — A — fc A + 



fe-i 



= fc 



1 f 1 
Y~77j J h k . x (u)&u- h k . x 



0+1 



fl - 0) 2 



24 «e(o,i) 



1 



= -fc(fc-l)(fc-2)A fc - 3 . 



Therefore, in order to have the inequality (2.15) hold, it is sufficient to make 
(fc- 2)A fc - 3 > 



fc-2 / ^ x fe-2 
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which is equivalent to 



A > 



k-2 



for all fc > 3. Since 



\ 



k-2 



k-2 /^\ fc-2 



< 



yfc — 2 V 2 



fc-2 



as fc — > oo, it follows that when A > | the inequality (2.15) holds. This means that 
when A > | the negative of the function (2.11) is completely monotonic on (0,oo). 
The proof of (2) in Theorem 1 is complete. 

3. Applications 

In this section, we apply the complete monotonicity of ±H\(x) to establish 
inequalities for bounding the gamma function, the n-th harmonic number X^fe=i \i 
and the factorial n\. 

Theorem 2. For x e (0, oo), the gamma function T(x + 1) can be bounded by 



(^) 



x + l/2 



exp 



1 



24(x + 1/2) 



< T(x + 1) 



< 



M^r 2 <^-^- 1 >)) ^ 



and 



x + 1/2 



x+l/2 



eXP >24 



2x 



12(ln7r - 1) 

x+l/2 



3(x + 3/2) 

/—z^ + i^v^^ ( 



< r(.x + 1) 
l 

+ 3/2) 



(3.2) 
(3.3) 



Proof. By (2.4), it is easy to see that 

lim H x {x) = 0. 

x^oo 

Moreover, it is immediate that 

lim H x (x) = — ( i + 12-121n7T 

x ^o+ 24 \ A 

By Theorem 1, it readily follows that when and only when < A < \ the function 
H\{x) is decreasing on (0, oo). So, we have 

< H x {x) < J-( I + 12- 121n7r 



24 V A 



that is, 



.x + l/2 



x+l/2 



1 



< V2~ir 



exp 
x + l/2 



24(x + A) 

x+l/2 

exp 



<r(x + i) 

1 + 12- 121n7r ] — 

24 I A x + A 



(3.4) 
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for < A < \ and x € (0,oo). The inequality (3.1) is proved. 

Similarly, when A > |, the function H\{x) is increasing on (0, oo), and the 
inequality (3.4) is reversed on (0, oo) for A > |. The inequality (3.2) follows. □ 

Theorem 3. For neN, the n-th harmonic number J2k=i J can ^ e bounded by 
ln ( n + I) + 24(n + iJW 

<ln ( n+ ^) + 24(n + 1 l/2) 2 + ^ (3 ' 5) 



ln ( n+ ^) + 24(n + 3/2) 2 +7< g^ 

<ln|n+i)H ; — 1 — tt + 1 - In - - — . (3.6) 

" V 2/ 24(n + 3/2) 2 2 90' 1 ' 

where 7 = 0.577 • • • stands for Euler-Mascheroni's constant. 

Proof. By Theorem 1 and the definition of completely monotonic functions, it fol- 
lows that 

(1) when < A < i, the function H' x (x) is increasing on (0, 00), 

(2) when A > |, the function H' x (x) is decreasing on (0,oo). 
Since 

and lim-r^oo H' x (x) — 0, it follows readily that 

ln ( X + + 24(x + 1/2)2 + 1 - 7 - ln I - k ~ ^ + 1} 



<ln(x+i)H ; — - — r (3.7) 

[ 2J 24(x + 1/2) 2 1 ' 



and 



l »{ X + l) + 24(^ + 3/2)2 <^ + 1 ) 

^ lll ( X+ + 24(, + 1 3/2)2 +1 - 7 - 1 4-^ (3 ' 8) 
for x e [1, 00). Taking x — n and using 

^(„+i) = £;I- 7) (3.9) 
fc=l fc 

see [1, p. 258, 6.3.2], in (3.7) and (3.8) give inequalities (3.5) and (3.6). □ 

We recall from [2, 22] that a function / is said to be logarithmically completely 
monotonic on an interval J C R if it has derivatives of all orders on / and its 
logarithm ln / satisfies 

(-l) fe [ln.f(x)]« > (3.10) 
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for k € N on /. 

Theorem 4. For x e (0, oo) and \>0, let 

e x T(x + l) 



{x + l/2)*+ 1 / 2 Cxp 24(x + X)' 
Then the function G\(x) has the following properties: 

(1) if and only if < A < \, the function G\(x) is logarithmically completely 
monotonia on (0,oo); 

(2) if the inequality (2.1) is valid, the reciprocal of the function G\(x) is loga- 
rithmically completely monotonic on (0,oo); 

(3) if X > |, a special case of the inequality (2.1), the reciprocal of the function 
G\(x) is logarithmically completely monotonic on (0,oo). 

Proof. This follows from the obvious fact that 

and the definition of logarithmically completely monotonic functions. 
Theorem 5. For neff, the factorial n! can be bounded by 



(3.11) 



□ 



n+1/2 X "+1/2 



exp 



1 



< V2^ 
and 



n + l/2\ n+1/2 



cxp 



24(n+ 1/2) 
1 



< n\ 



24 



12 3-ln7r + ln- 



27) 3(n+l/2) 



(3.12) 



n + 1/2 



-1/2 



exp 



< V2H 



121 3-ln7r + ln — 
27 

n + 1/2 N 



1 



exp 



5(n + 3/2) 
1 



< n\ 



24(n + 3/2) 



(3.13) 



Proof. Combining (3.3) and 



24 



— !— + 36-12 ln(27r) - 36 In - 
A + 1 2 



with Theorem 1 and the proof of Theorem 2 reveals 



x + 1/2 



x+1/2 



exp - 



1 



x + 1/2 



x+1/2 



24(x + 1/2) 
1 



exp 



24 



< T(x + 1) 
4 



12 3-ln7r + ln- 



27 J 3(:r+- 1/2) 



(3.14) 



and 



— fx + l/2\ x+1/2 (1 

T(x + 1) < >/2tt^ 



12 3 - In 7T + In — 
27 

x + l/2\ x+1/2 



+ 



1 



cxp - 



5(.t + 3/2) 
1 



< 



24(^ + 3/2), 



(3.15) 
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Letting x — n and using T(n+1) = n\ in (3.14) and (3.15) leads to inequalities (3.12) 
and (3.13). The proof is complete. □ 

4. Remarks 

In this section, we would like to comment some results above-presented. 

Remark 1. The inequality (2.14) and the inequality (2.15) for A > | may be 
rearranged as the double inequality 

1 x 2 x I x 2 



e x/2 24e^ x ~ e x - 1 ~ e x / 2 24e ax 

on (0, oo), where a = § and /3 = 5. This improves the right hand side and partially 
improves the left hand side of the double inequality 

4 < ~1T~~\ < 4n> x > ( 4 - 2 ) 

e x e x — 1 e x ' z 
in [23, p. 2550, Proposition 4.1]. 

We guess that the scalar a = | in (4.1) can be replaced by a smaller number, 
for example, 1, but the constant (5 = \ in (4.1) is the best possible. 

In [19], some related inequalities for the exponential function e x were constructed. 

In a subsequent paper, we will refine the right hand side of the double inequal- 
ity (4.1) and employ it to strengthen double inequalities for bounding Mathieu's 
series 

^ Or, 

%) = Efe?T7¥' r >° ( 4 - 3 ) 

n=l ^ ' 

and the like. For more information on bounding Mathieu type series, please refer 
to [16, 20, 23] and closely related references therein. 

Remark 2. There have been plenty of references devoted to bounding the n-th 
harmonic number Y^k=i j; for n G N, for example, [7, 10, 11, 24] and closely related 
references therein. 

Remark 3. Several inequalities for bounding the gamma function were also estab- 
lished and collected in [9, 12, 13]. See also [21, pp. 52-57] and lots of references 
cited therein. 

Remark 4. It was proved once again in [3, 8, 22] that the set of logarithmically 
completely monotonic functions is a subset of the completely monotonic functions. 
This implies that Theorem 4 is not trivial. 

Remark 5. In [14, 15], it was shown that the function 

. e x r(a; + l) „ „ 

9^) = (4-4) 

on the interval (max{0, — /3}, 00) for /3 G R is logarithmically completely monotonic 
if and only if j3 > 1 and that the function [ga^ix)]' 1 is logarithmically completely 
monotonic if and only if /3 < |. See also [21, pp. 53-54, Section 5.6]. Motivated by 
this and Theorem 4, we would like to ask a question: How about the logarithmically 
complete monotonicity of the function 

o 1 ^ e x r(a; + l) 1 u 
G\ u( x ) = 7 N , exp — — (4.5) 

on the interval (max{0, —A, — /i}, 00)? where A and /x are given real numbers. 
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